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We derive exact thermodynamic identities relating the average number of condensed atoms and 
the root-mean-square fluctuations determined in different statistical ensembles for the weakly inter- 
acting Bose gas confined in a box. This is achieved by introducing the concept of auxiliary partition 
functions for model Hamiltonians that do conserve the total number of particles. Exploiting such 
thermodynamic identities, we provide the first, completely analytical prediction of the microcanoni- 
Oh! cal particle number fluctuations in the weakly interacting Bose gas. Such fluctuations, as a function 

of the volume V of the box are found to behave normally, at variance with the anomalous scaling 
behavior V'^^^ of the fluctuations in the ideal Bose gas. 



00 



> 



0^ 
0^ 



After the experimental achievement of Bose-Einstein condensation in trapped, cold gases of weakly interacting 
neutral Bose atoms [0- a flow of experimental data on the statistical properties of the condensate fraction has 
become available for theoretical study Among others, the issue of the fluctuations in the number of condensed 
I atoms is of central importance, and it is foreseeable that it will become experimentally testable in the near future. 
' Experimentally, the trapped atoms are in a situation of almost complete isolation with respect to the outer envi- 
ronment surrounding the trap: the total number of particles is certainly conserved, while energy exchange with the 
0^ ' environment is small, in some cases even negligible. Therefore, the relevant statistical ensemble to consider when 
] attempting to compute the fluctuations should be the canonical or, even more appropriate, the microcanonical one. 
For the ideal Bose gas, the fluctuations of the condensate fraction were calculated in a series of papers Q- [ pd] , 
both in the microcanonical and the canonical ensembles by a variety of approximate methods: both numerical (based 
I on recurrence formulae and on contour integration) and analytic (asymptotic formulae based on the notion of the 

a Maxwell Demon ensemble), for various trapping potentials and in different number of dimensions. 
I , On the other hand, the issue of fluctuations in the weakly interacting Bose gas is not as well defined, both because 
'"^ ■ there is no unambigous way of defining the condensate fraction and because the approximate dynamics and energy 
^ ] spectra can be defined in a number of different ways, yielding different model-dependent predictions. 
O ■ Of the few existing attempts to address the problem, we should mention in particular the calculation of the 
I fluctuations of the interacting Bose gas trapped in a box, that was developed analytically in the canonical ensemble 
assuming the Bogoliubov number nonconserving description Jl^ . In this framework, it was predicted an anomalous 
scaling of the fluctuations with the volume V of the box, < 5~No >cn^ V^l"^ analogous to that of the ideal Bose gas. 
However, this result can be questioned on the grounds that it has been derived assuming the condensate to act as a 
particle reservoir, and moreover that in the actual experimental situations the total number of particles A'^ is a good 
quantum number. 

It is therefore sensible to set up model calculations that avoid these features, and to compare the results thus 
obtained. As a first step, the fluctuations were calculated analytically in the canonical ensemble and numerically 
in the microcanonical ensemble, assuming different standard approximate spectra (respectively, the one obtained in 
Hartree-Fock theory, the two-gas model obtained by reduction of the Bogoliubov approximation, and the spectrum 
obtained in lowest order perturbation theory). It was shown that the fluctuations are very sensitive to the details 
of the approximation, resulting in large discrepancies between the different models jl^ . On the other hand, they all 
satisfy the criterion of conservation of the total number of particles, while the lowest order perturbative spectrum 
includes three- and four-body contributions to the dynamics. 

In this letter we address the issue of computing analytically the fluctuations of the condensate fraction in the 
microcanonical ensemble for a specific number-conserving model interaction. By introducing a so-called auxiliary 
canonical partition function, defined below, we prove by saddle point evaluation and with the help of standard 
thermodynamic Legendre transformations that for all model Hamiltonians that do conserve the total number of 
particles, the microcanonical fluctuations can be expressed in terms of the canonical moments exactly as in the case 
of the microcanonical fluctuations in the ideal Bose gas j|] . 

On the other hand, the canonical moments are modified in the presence of the interaction, and so is the final analytic 
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expression of the microcanonical fluctuations, which we are able to compute. At smafl but finite couphng (scattering 
length) the fluctuations scale linearly with the volume V of the box, showing no anomalous behavior. In all actual 
calculations we assume for the interacting system the energy spectrum obtained in lowest order perturbation theory 

The Bose-Einstein condensation manifests itself in the macroscopic occupation of one of the eigenmodes. In 
general, the eigenfunctions of the single particle density matrix depend on the total number of atoms N, on the 
trapping potential, on the interatomic interactions and on the statistical ensemble of the system. The exceptionally 
simple, and hence the best studied, is the case of a system confined in the box with periodic boundary conditions. In 
this case the single particle density matrix is translationally invariant and periodic, so its spectral decomposition is 
just the Fourier series, and the index j is the quantized momentum p = {2fnT / L){nx,ny,nz), where L is the size of 
the box. 

We are interested in the statistics of a weakly interacting Bose gas for temperatures below the condensation tem- 
perature. In this regime, the system may be viewed as being composed of two macroscopic subsystems, a condensate 
of A^o ~ 0{N) particles occupying the p = state, and an excited part of N^x = Ylp^o particles. For isolated 
systems, particle number conservation implies Nex — N — Nq. 

After splitting the system into its condensed and excited parts we will have to define the approximate dynamics. 
This needs not to be specified as long as we are concerned with general thermodynamic arguments, but only in the 
actual calculations. 

One commonly made assumption is to neglect the interatomic interaction within the excited subsystem. Thus, the 
states of the excited atoms are still single particle states in a box, since the condensate is uniform. Accordingly, the 
energy of the excited subsystem is given by 

p2 

To fully define the problem, we need to specify the relation between the total energy of the system E and E^x, 
subject to the requirements that the total number of atoms A^ be a good quantum number, and that higher order 
contributions beyond the two-body ones be included. We will hence adopt the spectrum obtained in lowest order 



perturabtion theory 1 14 1 in its simplified form by neglecting the terms proportional to the product of excited occupation 
numbers: 

E{N, N,x) = a{N'' + 2A^A^e. - iV^J + E,x , (2) 

where 

mV 

and a is the scattering length. 

The above, like the other mentioned approximate spectra, depends only on the total number of atoms A^ and on 
the total number of excited atoms N^x (or, equivalently, on A^ and Aq), and not on the detailed distribution among 
excited states. 

To calculate the canonical expressions for mean values and fluctuations and to relate them to the microcanonical 
ones in an exact thermodynamic setting, we may now proceed as follows: we define an auxiliary canonical partition 
function for exactly A" particles, by introducing the chemical potential /iex pertaining to the subsystem of excited 
atoms: 

Zaux{N, (3, fiex) = TrN [cxp -(/3i/ - ^iexNex)] = ^ exp [~^iexNex]Zex{N, Nex,f3) ■ (4) 

The crucial observation is now that differentiating the free energy with respect to fj,ex and evaluating the derivative 
at flex = yields the canonical average < A'o >cn- Furthermore, since N^x commutes with the total Hamiltonian, 
second derivative evaluated at fi^x — yields the canonical fluctuations < <5^Ao >cn- 

To proceed, one determines flrst In Z^x by the method of the most probable value. This is done expressing Z^x in 
terms of the associated generating function, i.e. the grand-canonical partition function 'E.^x- As explained below, this 
introduces a further chemical potential /i^, not to be confused with the chemical potential fi^x of the excited atoms 
introduced in Eq. (1). The saddle point value jig of this chemical potential is flxed by the value Nex of Nex for which 
In Zex acquires its maximum. Standard calculations lead to the following expression for the most probable value: 
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lnZe,(7V,7Ve„/3) = -pE{N,Ne^) + ji.N,, - V 



3/2 



35/2(exp-/Xs) 



(5) 



where (75/2 is the Bose function of order 5/2 and /Is — 2j3a{N — Nex) is positive defined. To simpUfy notations, we 
have absorbed a factor (3 in the definition, so that fls is adimensional, in contrast to the chemical potential ^ex- 
Saddle point evaluation of the total free energy at the equilibrium value N^x can now be performed to yield 



lnZauxiN,P,^^ex) = lnZe,(iV,iVe„/3) +/3Me.A^e. + ^ )^ ' 
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(6) 



We see by immediate inspection that the first derivative with respect to fiex evaluated at the equilibrium value 
Mea: = yiclds the canonical mean value < N^x >cn— Nex, while the second derivative (again evaluated at ^ex — 0) 
yields the canonical root-mean-square fluctuations < S'^N^x >cn— d'^lnZex/dN^x- 

The implicit equation satisfied by the canonical mean number of condensed atoms < Nq >cn= Nq is Nex — N — Nq, 
where the latter term reads: 
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where the Bose-Einstcin function gn{s) — X^i^i s'' /I". In the same approximation the canonical fluctuations read 



< S^No >cn = 
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^0 4sinh2 [f (£i + 2a{N - Nex)] 
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(8) 



These expressions, as it should be, coincide with the ones already obtained by Idziaszek et al. |T^] in the standard 
canonical setting. Now, by applying twice the Legendre transformation we can express the microcanonical thermody- 
namic derivatives (fundamental variables: E,N) in terms of the canonical ones (fundamental variables: P,N). One 
has: 
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The above expression for the microcanonical fluctuations can be recast by simple thermodynamic manipulations in 
a more transparent form. By noting that 
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we can write 



< S^No >r. 



< S^No >c 



[< SNqSE >enY 

< S^E >™, 



(11) 



Relation (11) is our first main result. It coincides with the formula derived by Navez et al. Q for the ideal Bose 
gas. We thus find that the general structure of thermodynamics and thermodynamical fluctuations is not modified by 
the presence of a model interaction that conserves the total number of particles. In other words, our analysis satisfies 
the general principle that the thermodynamic relations between different statistical ensembles should be independent 
of the dynamics (the Hamiltonian) , a principle that is violated when resorting to Bogoliubov theory. 
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Direct analytic evaluation of the microcanonical fluctuations is now possible by computing the derivatives of the 
approximate free energy (Eq. (6)) appearing in Eq. (9). The final result is: 

= -IM^. (12) 

where 
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3ff3/2(exp i~^ls)) + 4Q/3A^ogi/2(exp (-^^)) 
2/3[l-2A(/?)a/?<7i/2(exp (-/i,))] 



G(/3,/2,) = ^A(/3) 
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/32 a/3 



53/2(exp {-[Is)) + + ^^(/3)35/2(exp (-/^s)) ■ (13) 



The above expression is our second main result. It generalizes the microcanonical fluctuations of the number of 
condensed atoms in the ideal Bose gas to the interacting case. For small values of the coupling a (i.e. of the scattering 
length a) we can exploit the Robinson-Kac representation of the Bose g functions up to terms linear in {—(i^^l'^. 
In this situation, the canonical fluctuations behave normally (scaling linearly with V) as already shown by Idziaszek 

et al. while the second term in the right-hand-side of Eq. (12), up to order ^lJ"^ reduces to 



3A(/3) 
535/2(1) 



53/2(1) + 2V^(/2,)i/2]' . (14) 



It is straightforward to verify that this term too scales linearly with V ^ also in the thermodynamic limit — > 00, 
iV — > 00, N jV = const. Therefore the total microcanonical fluctuations scale linearly with the volume of the box, 
and show no anomalous behavior for any finite value of the coupling. Obviously, in our model the thermodynamic 
limit does not commute with the weak coupling limit, and to recover the anomalous scaling of the fluctuations in the 
ideal Bose gas one must first perform the thermodynamic limit and then go to the limit a — » of vanishing scattering 
length. 

In Fig. 1 the microcanonical root-mean-square fluctuations Eq. (14) are plotted (dotted line) and compared to the 
canonical results Eq. (14) of Ref. |0| (dot-dashed line), to the canonical fluctuations computed in the framework of 
Bogoliubov theory, Eq. (8) of Ref. [Il2| (dashed line) and to the ideal Bose gas behavior (solid line) . The fluctuations 
are plotted as a function of the temperature for N = 10"* particles and a/ L = 5 x 10""*. The temperature is measured 
in units given by the spacing between the two lowest levels in the 3D box: A = {2Tifif' / {2raL'^) . The inset displays 
results for N = 10® and a/L = 5 x 10~^; the curves correspond to our microcanonical result (dashed), the ideal Bose 
gas (solid), and the canonical results of Ref. |l2) (dashed). The canonical results of Ref. jl^ are not shown in the 
inset because with increasing number of particles they become indistinguishable from the microcanonical prediction, 
as expected. 

We see that our result (14) is in marked disagreement with the predictions obtained in [ p^ . The emerging picture 
of ideal versus weakly interacting Bose systems is then the following. In the framework of Bogoliubov approximation: 
the anomalous scaling behavior of the ideal gas fluctuations is preserved; the general thermodynamic relation between 
ideal gas fluctuations computed in different ensembles is lost. In the framework of number conserving approximations: 
the scaling behavior of the ideal Bose gas fluctuations is modified; the ideal Bose gas thermodynamic relations are 
preserved. 

To summarize, we have performed a model calculation of the microcanonical fluctuations of the weakly interacting 
Bose-Einstein condensate, allowing for analytical solution, to provide a first prediction that may be tested against 
experiment. We stress this point, since in the real experimental situations the trapped gases are usually in conditions 
of fantastic isolation, and thus the microcanonical result should be the one to be compared with actual experimental 
observations (when they will become available). 



At variance with the canonical fluctuations predicted in the framework of Bogoliubov theory |12|, we find no 
anomalous scaling of the microcanonical fluctuations for any finite value of the scattering length. On the other hand, 
exact knowledge of the full energy spectrum is available only for the one-dimensional chain of interacting bosons . 
In this case, comparison of the exact solution with Bogoliubov theory shows that the Bogoliubov spectrum accounts 



4 



only for some subset of the low-lying elementary excitations. The situation for Bogoliubov theory is likely to be 
even worse in higher dimensions. Furthermore, phonons are collective excitations which do not necessarily change 
the number of particles in the condensate. It is therefore hardly conceivable that they should play an important role 
in determining the equilibrium fluctuations in the occupation number, and may therefore be safely ignored in first 
approximation. 

Some of the most subtle aspects of the dynamics and statistical mechanics of many-body interacting Bose systems 
appear in the above discussion: care is needed in handling the different limits, whose order cannot in general be 
interchanged, while the non-trivial thermodynamic identities, relating the different statistical ensembles at the level 
of fluctuations, are found to hold true just as in the case of the ideal Bose gas 

Looking at the future perspectives, we expect that achieving a deeper understanding of the thermodynamic prop- 
erties in exactly solvable instances such as the Lieb model in one dimension would be of great help, and studies 
addressing this problem are currently under way. On the experimental side, we hope that the measurement of the 
second-order correlation function (?2(?'i7?'2) will soon become available, allowing for a direct test of the different 
theoretical predictions. 
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FIG. 1. Root-mean-square fluctuations of the ground state occupation number for N — 10'* and a/L = 5 x 10~^. 
Displayed are (i) our microcanonical result Eq. (14) (dotted), (ii) Eq. (14) of Ref. |l3) (dot-dashed), (iii) the ideal 
Bose gas (solid line), and (iv) Eq. (8) of Ref. ||T^ (dashed). Inset: root-mean-square fluctuations for N = 10^ and 
a/L = 5 X 10~^. Displayed are (i) our microcanonical result (dotted), (ii) the ideal Bose gas (solid line), and (iii) Eq. 
(8) of Ref. H (dashed). 
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